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ABSTRACT 

We present high resolution (i?=80,000) spectroscopy of [Nell] emission from 
two young stars, GM Aur and AA Tau, which have moderate to high inclinations. 
The emission from both sources appears centered near the stellar velocity and 
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is broader than the [Nell] emission measured previously for the face-on disk 
system TW Hya. These properties are consistent with a disk origin for the 
[Nell] emission we detect, with disk rotation (rather than photoevaporation or 
turbulence in a hot disk atmosphere) playing the dominant role in the origin 
of the line width. In the non-face-on systems, the [Nell] emission is narrower 
than the CO fundamental emission from the same sources. If the widths of both 
diagnostics are dominated by Keplerian rotation, this suggests that the [Nell] 
emission arises from larger disk radii on average than does the CO emission. The 
equivalent width of the [Nell] emission we detect is less than that of the spectrally 
unresolved [Nell] feature in the Spitzer spectra of the same sources. Variability 
in the [Nell] emission or the mid-infrared continuum, a spatially extended [Nell] 
component, or a very (spectrally) broad [Nell] component might account for the 
difference in the equivalent widths. 

Subject headings: (stars:) circumstellar matter - - (stars:) planetary systems: 
protoplanetary disks — stars: pre-main sequence — (stars: individual) AA Tau, 
GM Aur 



1. Introduction 

The [Nell] 12.8 /um emission line has been suggested as a potential new probe of the 
planet formation region of circumstellar disks. It is a potentially powerful diagnostic because 
the neon in disks is expected to be fully in the gas phase and in atomic form. In addition, the 
[Nell] 12.8 fj,m line probes warm, ionized gas, conditions which are believed to characterize the 
upper atmosphere of the inner disks surrounding classical T Tauri stars (Glassgold, Najita, & 
Igea 2007; Meijerink, Glassgold, & Najita 2008) and disk photoevaporative flows (Alexander 
2008). Because it is sensitive to low column densities of gas, [Nell] may also be a useful 
probe of residual gas surrounding weak-line T Tauri stars or gas in the optically thin regions 
of transitional disks. Glassgold et al. (2007) predicted that the inner regions (< 20 AU) 
of classical T Tauri disks that are irradiated by stellar X-rays would produce strong [Nell] 
emission that could be detected with the Spitzer Space Telescope. Comparably strong [Nell] 
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emission has indeed been detected, apparently commonly, in Spitzer spectroscopy of T Tauri 
stars (Pascucci et al. 2007; Lahuis et al. 2007; Ratzka et al. 2007; Espaillat et al. 2007). 

While the rough agreement between the predicted and observed [Nell] line strengths 
supports the interpretation of a disk origin for the emission, stronger confirmation can be 
obtained from a study of resolved [Nell] line profiles. Herczeg et al. (2007) previously re- 
ported high resolution spectroscopy of [Nell] emission from TW Hya. TW Hya is a face-on 
disk system, with an inclination of i — 4 degrees for the inner disk (Pontoppidan et al. 2008; 
see also Qi et al. 2004). It is also a "transition object", a T Tauri star whose spectral energy 
distribution (SED) indicates that the disk continuum is optically thin within a given radius. 
In the case of TW Hya, modeling of the SED indicates that the disk is optically thin in the 
continuum (A > 1 /im) within ~ 4 AU (Calvet et al. 2002; cf. Ratzka et al. 2007). 

In their study of the [Nell] emission from TW Hya, Herczeg et al. (2007) found that 
the emission was centered at the stellar velocity, consistent with a disk origin. However, 
the emission was also significantly broader (FWHM ~ 21kms~ 1 ) than has been predicted 
theoretically for a near face-on disk (FWHM < 5kms _1 ; Meijerink et al. 2008; Glassgold 
et al. 2007). Herczeg et al. suggested that the observed line width could be explained if 
the emission arises from (1) a rotating disk at much smaller disk radii (~ 0.1 AU) than in 
the model; (2) a disk with a significant turbulent velocity component; or (3) a face-on disk 
undergoing photoevaporation. While each of these scenarios might account for the observed 
emission from a face-on system such as TW Hya, Herczeg et al. (2007) proposed that these 
scenarios could be tested by observing [Nell] emission from higher inclination disk systems. 

A more recent paper (van Boekel et al. 2009) illustrates another possible origin for the 
[Nell] emission from T Tauri stars. High resolution spectroscopy of the [Nell] emission from 
the T Tau triplet reveals that the emission is spatially extended and associated with a known 
outflow in the system. This suggests that [Nell] emission may arise primarily in outflows 
rather than disks in systems with strong outflow activity. 

To investigate the scenarios proposed by Herczeg et al. (2007) for the origin of [Nell] 
emission in T Tauri stars, we present here high resolution spectroscopy of the [Nell] emission 
from AA Tau and GM Aur, T Tauri stars that are not known to show strong outflow activity. 
Apart from its high inclination (i = 75 degrees; Bouvier et al. 1999), AA Tau is a typical 
classical T Tauri star. GM Aur is an actively accreting transition object that is viewed at an 
intermediate inclination (i = 54 degrees; Simon et al. 2000). The SED of GM Aur indicates 
that the disk is optically thin in the mid-infrared continuum, suggesting that the disk is 
devoid of small dust grains in the radial range 5-24 AU (Calvet et al. 2005). 
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2. Observations and Data Reduction 

The observations were carried out using TEXES in its high-resolution mode (R = 
80, 000) with a 0.6" wide slit on the Gemini North telescope on 29-30 October 2007 un- 
der program ID GN-2007B-C-5. The background sky emission was removed by nodding 
the source along the slit and subtracting adjacent nod positions. We corrected for telluric 
absorption with observations of a CMa obtained at a similar airmass as that of the science 
targets (the difference in airmass was < 0.2). The spectra were flux calibrated using the 
low resolution Spitzer Space Telescope spectra reported by Furlan et al. (2006). Roughly 
every 10 minutes, we took a series of calibration frames including blank sky and an ambi- 
ent temperature blackbody. The blackbody observations were used for flatfielding the data, 
while the sky emission line observations were used for wavelength calibration. The standard 
TEXES pipeline (Lacy et al. 2002) which produces wavelength-calibrated one dimensional 
spectra, was used to reduce the data. The wavelength solutions have an accuracy in velocity 
of ~ 1 kms -1 . 

One challenge associated with these observations is that the [Nell] line at 12.81 //m falls 
in a gap between two TEXES spectral orders. By tilting a mirror within the instrument, we 
are able to cover the gap by shifting the spectrum on the detector and thereby accessing the 
[Nell] line at the long wavelength end of one order or the short wavelength end of the next 
order. In the first case, the blaze efficiency rises toward the blue side of the line, while in 
the second case, it rises toward the red side of the line. 

During the night of 29 October, we shifted the optics to enhance the sensitivity on the 
red side of the [Nell] line. We spent 2590 seconds of on-source integration time on GM Aur 
and 3238 seconds on AA Tau in this instrument setup. The observed GM Aur line appeared 
narrow, symmetric, and centered at the stellar velocity while the emission from AA Tau was 
broad and redward of the stellar velocity. This led us to search for a line profile component 
that might be present blueward of the AA Tau stellar velocity on the following night. We 
spent 2331 seconds of on-source integration time on AA Tau on the night of 30 October 
2007, with our optics shifted to enhance the sensitivity on the blue side of the [Nell] line. 

The AA Tau data from the two separate nights were combined by interpolating onto 
a common wavelength scale and, in the spectral regions where data from both nights were 
available, weighting the contribution from each night. The weighting was used to account 
for the varying noise across the spectra in the overlap region; the variation resulted primarily 
from the rising or falling blaze efficiency. We therefore weighted each data point inversely as 
the square of the noise, which was calculated assuming that photon noise dominates. This 
weighting is appropriate if the noise in the two spectra being combined add in quadrature. In 
practice, the exact value of the weighting factor had little effect on the combined spectrum. 
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3. Results 

Figures 1 and 2 show the resulting [Nell] spectra of AA Tau and GM Aur. The spectra 
shown have been smoothed by 3 pixels. Because 3 pixels is the width of a spectral reso- 
lution element at the wavelength of [Nell], the smoothing has little impact on our results. 
Although the larger noise on the blue side of the [Nell] line in AA Tau makes it difficult 
to characterize the line profile with certainty, the emission appears to be centered near the 
stellar heliocentric velocity (vhciio = 16.5 km s -1 ; Bouvier et at. 1999) and is broad, extending 
both redward and blueward of the line center with FWHM ~ 70kms _1 . The line appears 
approximately double-peaked, with a blueshifted peak possibly further from the stellar ve- 
locity and possibly weaker than the redshifted peak. The emission redward of the stellar 
velocity has an equivalent width of 6.6 ± 1.4 A. 

In comparison, the spectrally unresolved [Nell] feature detected in the i?=600 Spitzer 
IRS spectrum of AA Tau (Carr & Najita 2008) has an equivalent width of approximately 
24 A. If the emission blueward of the stellar velocity in the TEXES spectrum has an equiva- 
lent width equal to that of the redward emission, the total equivalent width (red and blue) 
is ~ 0.55 of the equivalent width seen in the Spitzer spectrum. Such a difference might 
arise if either the mid-infrared continuum or the [Nell] emission is time variable. Another 
possibility is that the Spitzer spectrum (4.7" slit width) includes spatially extended [Nell] 
emission that is excluded in the narrower slit (0.6") of the TEXES observation. Spatially 
extended [Nell] emission can arise in outflow sources (e.g., Neufeld et al. 2006; van Boekel 
et al. 2009). 

Spectrally unresolved emission from other lines might also enhance the equivalent width 
of the [Nell] feature in the Spitzer spectrum. However, known emission features are unlikely 
to contribute at this wavelength. The synthetic spectrum of AA Tau presented in Carr & 
Najita (2008) provides a good fit to the molecular emission detected in the Spitzer spectrum 
of that source. When examined at high spectral resolution, the synthetic spectrum reveals 
no strong lines of H 2 0, HCN, or C2H2 within 200 km s -1 of the [Nell] line. We therefore 
interpret the emission features detected with TEXES as [Nell]. 

The [Nell] emission from GM Aur has a component of modest width (FWHM ~ 
14kms _1 ) that is centered near the systemic velocity determined from measurements of 
the surrounding molecular disk (vhciio — 14.8 km s -1 ; Dutrey et al. 1998; Simon et al. 2000). 
There is possible evidence for an additional emission component centered ~ 40kms _1 red- 
ward of the stellar velocity, although it is difficult to be certain because of the limited 
signal-to-noise ratio of the spectrum. The 40kms~ 1 component might arise from either in- 
falling gas located close to the star or gas in rotation about the star. In the latter case, we 
would expect to see an additional emission component with a similar ~ 40kms _1 velocity 
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offset on the blue side of the line. Because we observed the [Nell] emission in only one order, 
the GM Aur spectrum has much higher noise on the blue side of the line (see §2), and it is 
impossible to determine whether a corresponding blueshifted component is present. A more 
complete study of the [Nell] emission from GM Aur is needed to distinguish between these 
possibilities. 

The emission feature centered near the stellar velocity has an equivalent width of 18±2A, 
and the feature ~ 40kms~ 1 redward of the stellar velocity has an equivalent width of 8±2A. 
In comparison, the equivalent width of the spectrally unresolved [Nell] feature detected in 
an i?=600 Spitzer IRS spectrum of GM Aur (Najita et al., in preparation) is approximately 
70A, almost three times the combined equivalent width of the two [Nell] emission features 
detected at high spectral resolution. As discussed above, variability in the [Nell] line or the 
mid-infrared continuum, alternatively spatially extended [Nell] emission, might account for 
the difference in the equivalent widths. We are also insensitive to very broad [Nell] emission, 
particulary given the limited signal-to-noise ratio in the continuum. A line 10% above the 
marked continuum in Figure 2 and 300 km s -1 wide, which is within a Spitzer resolution 
element, would have an equivalent width of 13 A, similar to the equivalent width reported 
for the [Nell] line. Such a broad component would be difficult to detect with our data. 

Another possibility is spectrally unresolved emission from other lines. In the Spitzer 
spectrum of GM Aur, the emission feature detected at the wavelength of [Nell] is broader 
than an unresolved line. Perhaps the feature includes a contribution from lines other than 
[Nell]. Mid-infrared molecular features such as H2O, HCN, or C2H2 would contribute very 
negligibly to our TEXES observations, because any such emission features, if present in the 
Spitzer spectrum of GM Aur, are much weaker than in AA Tau. We therefore interpret the 
emission features detected with TEXES as [Nell] for the purpose of this paper. 

Thus, as in the case of TW Hya (Herczeg et al. 2007), the [Nell] line profiles of AA Tau 
and GM Aur are consistent with emission centered near the stellar velocity. In comparison 
with the intrinsic (deconvolved) line width of 21kms _1 (FWHM) obtained by Herczeg et al. 
(2007) for TW Hya, the width of the [Nell] emission from AA Tau is significantly broader 
(~ 70kms _1 ). While the [Nell] emission from GM Aur has a velocity component centered 
on the stellar velocity that is more similar in width (~ 14kms _1 ) to the [Nell] emission 
from TW Hya, the [Nell] emission may also include a redshifted emission component that 
extends to > 50kms _1 of the stellar velocity. 
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4. Discussion 
4.1. Comparison to Profiles of Other Emission Lines 

4-1.1. Ha Profiles 

The double-peaked shape of the [Nell] line profile of AA Tau is similar to the shape 
of the Ha line profile of AA Tau (Bouvier et al. 2007, 2003, 1999; Alencar & Basri 2000; 
Edwards et al. 1994). In T Tauri stars, the Balmer lines are believed to form in stellar 
magnetospheres (e.g., Calvet & Hartmann 1992; Edwards et al. 1994; Muzerolle et al. 1998). 
Although both the Ha and [Nell] profiles of AA Tau are double-peaked, the Ha profile is 
also much broader (FWHM ~ 280 km s" 1 ; Bouvier et al. 1999). 

For GM Aur, the Ha emission is centrally peaked (Edwards et al. 1994), as is the 
central component of its [Nell] emission. However, the Ha emission (FWHM of ~ 220 km s -1 ; 
Edwards et al. 1994) is also much broader than the [Nell] emission (~ 14kms _1 FWHM); the 
half-width at zero intensity (HWZI) of the Ha line (~ 500 km s -1 ) is also much broader than 
the velocity extent of the redward emission of the [Nell] emission (~ 50kms _1 ). While the 
Ha and [Nell] profiles have some morphological similarities, the difference in the velocity 
widths suggests that the [Nell] emission arises primarily from a different region than the 
stellar magnetosphere. 

4-1.2. CO Fundamental and UV H 2 Profiles 

The CO fundamental lines from AA Tau are also double-peaked, similar to the shape of 
the [Nell] line in AA Tau. CO fundamental emission is believed to probe the inner regions 
of T Tauri disks (see Najita et al. 2007 for a review). Figure 3 compares the [Nell] profile of 
AA Tau with the average line profile of CO emission from AA Tau that was obtained by J. 
Carr, J. Najita, and N. Crockett on the night of 26 Nov 2004 using NIRSPEC on the Keck 
II telescope on Mauna Kea. The line profile shown is the average of 4 i>=l-0 lines in the 
4.9 fxm region. Regions of poor telluric correction (< 85% transmission, as determined from 
observations of the telluric standard) were excluded from the average. The double-peaked 
CO line profile has a FWHM of ~ 140 km s -1 and a peak-to-peak separation of ~ 80kms _1 . 
Thus, the CO profile is broader than the [Nell] line profile, but narrower than the Balmer 
line profile. A preliminary analysis of the basic properties of the CO fundamental emission 
from AA Tau finds that the emission is optically thick and has an average temperature of 
~ 900 K (Carr & Najita 2008). A more detailed analysis of the CO emission will be presented 
in a future publication. 
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CO fundamental emission has been reported from GM Aur by Salyk et al. (2007), al- 
though it was not detected in an earlier study (Najita et al. 2003). Because of the low infrared 
continuum excess of GM Aur compared to that of AA Tau (GM Aur is a transition object), 
the stellar photosphere makes a more significant contribution to the GM Aur spectrum than 
it does in the case of AA Tau. We therefore first corrected for stellar photospheric absorption 
in the Salyk et al. (2007) spectrum before constructing the average CO line profile. Following 
the approach described in Najita et al. (2008), we created a synthetic stellar spectrum that 
is appropriate for the stellar spectral type of the source (K3V; Herbig & Bell 1988). We also 
used a stellar rotational velocity of vsini = Skms^ 1 (cf. Hartmann et al. 1986). We found 
a good fit to the observed structure in the continuum (i.e., in the region away from the 1-0 
CO lines) with a veiling of 1 at 4.7 /im. This value for the veiling at 4.7 ^m is consistent with 
that implied by the SED of GM Aur (Furlan et al. 2006). We subtracted the veiled stellar 
continuum from spectrum, and we then added back an equivalent (featureless) continuum 
in order to show the strength of the emission relative to the continuum. 

Figure 4 compares the [Nell] profile of GM Aur with the average CO line profile con- 
structed from the resulting spectrum in the region of the v =1-0 P8 through P12 transitions. 
Regions of poor telluric correction (< 75% transmission) were excluded from the average. 
The average CO line profile differs from that reported by Salyk et al. (2007) in that there 
is no significant central dip in the average line profile, because the stellar photospheric con- 
tribution to the line profile has been removed. Because the average line profile is now more 
centrally peaked, the FWHM of the average line profile (~ SOkms^ 1 ) is also smaller than 
the value of ~ 50kms _1 reported by Salyk et al. (2007). 

The [Nell] and CO line properties can also be compared with those of UV fluorescent 
H 2 emission, another line diagnostic that is believed to probe the conditions in inner cir- 
cumstellar disks (e.g., Herczeg et al. 2006 and references therein; see e.g., Najita et al. 2007 
for a review). In the case of TW Hya, Herczeg et al. (2007) found that the [Nell] emission 
width (21kms _1 ) is broader than the UV H 2 emission (14kms _1 ), which is broader than 
the CO fundamental emission (~ Skms^ 1 ; Salyk et al. 2007). A similar decreasing sequence 
is inferred for the nominal temperatures of the emitting gas that gives rise to these features: 
~ 4000 K for [Nell] (Glassgold et al. 2007); ~ 2000 K for H 2 (Herczeg et al. 2004); and 
~ 800 K for CO (Salyk et al. 2007). Such a progression might occur if disk atmospheres 
experience transsonic turbulence, which is enhanced in the warmer, atomic regions higher 
up in the disk atmosphere. While this interpretation might be considered for a face-on sys- 
tem such as TW Hya, the progression of line widths in TW Hya, from broader [Nell] lines 
to narrower H 2 and CO lines, would not be expected to hold in higher inclination systems 
where disk rotation would play a larger role. 
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Indeed, for AA Tau, we find that the width of the CO emission is instead broader than 
the [Nell] emission. For GM Aur, the CO line (FWHM ~ SOkms -1 ) is also broader than 
the central component of the [Nell] line (FWHM ~ 14kms _1 ); the HWZI of the average CO 
line profile could be as large as ~ 50kms -1 (Salyk et al. 2007), comparable to the velocity 
extent of the red component of the [Nell] emission (~ 50kms -1 ). Thus, the CO emission 
is broader or comparable to the width of the [Nell] emission in the two sources studied. If 
both diagnostics are dominated by disk rotation, this suggests that the [Nell] emission arises 
from larger disk radii on average than the CO. 

This is consistent, in general, with models of the ionization and thermal structure of 
T Tauri disks (e.g., Glassgold et al. 2004, 2007; Meijerink et al. 2008). In the Glassgold 
et al. models, the surface of the disk is heated and ionized by stellar X-rays, producing a 
low column density (N H ~ 10 20 cm~ 2 ) surface layer of hot (~ 4000 K) gas that extends to 
> 10 AU. Significant [Nell] is found to emerge from the disk surface region within ~ 20 AU 
(Glassgold et al. 2007; Meijerink et al. 2008). In contrast, CO becomes abundant at much 
larger disk vertical column densities (Ah > 10 20 cm -2 ), and the warm temperatures and high 
densities needed to produce the CO emission (> 500 K) likely limits the emission to smaller 
disk radii, within a few AU. 

Although for TW Hya the width of the UV H_2 emission is broader than the CO fun- 
damental emission, the inverse is true for other sources with UV H2 and CO fundamental 
emission line widths that have have been reported in the literature (see Table 1). Thus, as 
might be expected, the progression of line widths found for TW Hya among the [Nell], H 2 , 
and CO emission lines does not appear to be typical of T Tauri stars. 

4.2. Origin of [Nell] Emission 

Herczeg et al. (2007) offered three possible interpretations for the origin of the width 
of the [Nell] line observed from TW Hya: (1) Keplerian rotation from the very inner disk 
region, (2) transsonic turbulence in a ~ 10, 000 K disk atmosphere, or (3) photoevaporation 
at ~ lOkms -1 from both faces of a face-on gaseous disk that is optically thin in the mid- 
infrared continuum. While each of these explanations could, in principle, account for the 
observed ~ 21kms _1 line width of TW Hya given its low inclination, they would produce 
different line profiles for systems viewed at higher inclinations. 

At higher inclination, a line profile dominated by disk rotation would remain symmetric 
about the stellar velocity but would become broader. In comparison, a line profile dominated 
by turbulence would show little difference when viewed at higher inclination. In the case of 
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a photoevaporative flow from a disk that is optically thick in the mid-infrared continuum, 
the disk would occult the receding flow, producing a line profile that is blueshifted by up to 
~ 10 kins" 1 , the nominal velocity of a photoevaporative flow (e.g., Font et al. 2004). 

A physical picture analogous to this last scenario is invoked to account for the high 
velocity blueshifted components of the [01] 6300A line emission from T Tauri stars (e.g., 
Edwards, Ray, & Mundt 1993); although the high velocities observed in the [OI] case indicate 
a dynamical (rather than thermal) ejection process. Font et al. (2004) studied whether a 
lower velocity thermally-driven disk photoevaporative flow could account for the properties 
of the low velocity component of the [OI] emission. A more recent paper explores the disk 
photoevaporative flow hypothesis specifically for [Nell], in both the case of a continuous disk 
and a disk with an inner hole (Alexander 2008). In a system viewed edge-on, the [Nell] line 
profile is predicted to be dominated by Keplerian disk rotation and to be centered on the 
stellar velocity. At lower inclinations, the [Nell] line profile is predicted to be blueshifted 
from the stellar velocity by 5 — lOkms" 1 . 

Of the three explanations invoked to explain the width of the [Nell] line of TW Hya, 
only disk rotation could plausibly account for the broad profile of the [Nell] emission from 
AA Tau. The width of the emission greatly exceeds the width expected from either disk 
turbulence or photoevaporation. The apparent double-peaked line profile might indicate 
that the emission arises from a limited range of disk radii, similar to the interpretation given 
to the double-peaked CO overtone lines in sources such as WL16 (Carr et al. 1993; Najita 
et al. 1996). For Keplerian rotation about a 0.7 M Q star viewed at an inclination of 75 
degrees, the redward velocity extent of the [Nell] emission and the redward emission peak 
(to +45kms _1 and at ~ 18kms _1 from the stellar velocity, respectively) would indicate 
emission extending from an inner radius of ~ 0.3 AU out to an outer radius of ~ 1.8 AU. 

The narrow (~ 14kms" 1 FWHM) component of the [Nell] emission from GM Aur that 
is centered at the stellar velocity could plausibly be explained by turbulent broadening. It 
might also arise from a rotating Keplerian disk at radii beyond several AU given the stellar 
mass (0.84 — 1.0M Q ) and inclination (i = 54 degrees) measured for the system (Simon et 
al. 2000; Dutrey et al. 2008). The lack of a velocity shift of this component from the stellar 
velocity argues against photoevaporation playing a significant role in the emission (cf. Font 
et al. 2004; Alexander 2008). The high velocity of the redshifted component could plausibly 
be accounted for by disk rotation if an (undetected) blueshifted component is also present. 
The ~ 40kms _1 velocity shift measured for this component would correspond to gas in 
Keplerian rotation at an orbital radius of 0.3 AU. 

In comparison, Glassgold et al. (2007) predicted that an X-ray irradiated disk atmo- 
sphere would produce [Nell] emission from radii within ~ 20 AU. [Nell] emission is expected 
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to arise from a similar range of disk radii in the case of a disk irradiated by stellar UV 
photons (Gorti & Hollenbach 2008; Alexander 2008). In an updated study that extends the 
Glassgold et al. (2007) calculation inward to 0.25 AU, Meijerink et al. (2008) showed the 
predicted [Nell] emissivity per unit radius, P, on a velocity scale i>k, where vk is the Keple- 
rian rotational velocity at a given radius. This radial contribution function corresponds to 
a [Nell] line profile with a FWHM approximately equal to 0.75 v\ sin i, where v\ is the Ke- 
plerian velocity at 1 AU. The line width is smaller than the width of the P(v K ) distribution 
because the line profile accounts for how the flux from a given annulus, when viewed at an 
inclination i, is distributed over a range of velocities between and vk- 

For the stellar mass and inclination of GM Aur, the width of the [Nell] line is expected to 
be 17kms _1 , similar to the observed 14kms _1 width of the central [Nell] emission component 
of GM Aur. Figure 5 compares the observed [Nell] profile of GM Aur with a model profile 
constructed from the P(vk) distribution of Meijerink et al. (2008), where the emission extends 
from an outer radius of ~ 25 AU, beyond which the [Nell] emission is insignificant, in to 
the minimum radius of 0.25 AU studied by those authors (short dashed blue line). The 
red-shifted [Nell] emission component extends to higher velocities, ~ 40kms~ 1 , and is not 
well accounted for by the model profile. For the i = 75 degree inclination of AA Tau, the 
Meijerink et al. (2008) results would suggest a width for the [Nell] line of 15 — 18kms _1 for 
a stellar mass of 0.5 — O.7M . The observed width of the [Nell] emission from AA Tau is 
significantly broader than this value. 

The difference between the observed profiles and the theoretical predictions may be due 
to the simple physical picture adopted in the disk emission models; real systems are likely 
to be more complex. Glassgold et al. (2004) adopted the conventional (flared, azimuthally 
symmetric) disk geometry of D'Alessio et al. (1999). In contrast, in AA Tau the modulation 
of its photometric and spectroscopic properties suggest that the inner disk edge (near the 
corotation radius) has a large scale height over a range in azimuth (Bouvier et al. 1999, 2003, 
2007). Periodic variations in the optical light curve occur in phase with enhanced redshifted 
absorption in the Balmer line profiles. These and other observed properties have been 
interpreted as evidence for an optically thick occulting screen produced by a magnetically- 
warped dusty inner disk edge. In a system with an inflated inner disk edge, more of the 
stellar ionizing photons (X-rays or UV) may be deposited close to the star, resulting in 
[Nell] emission from smaller radii and higher disk rotational velocities. 

GM Aur is also an unusual source in that, like TW Hya, it is a transition object. Based 
on the shape of the SED, Calvet et al. (2005) estimate that the disk is devoid of dust in 
the radial range 5-24 AU, with a small amount of dust present within 5AU. Earlier studies 
of the SED found a smaller optically thin region out to ~ 4AU (Rice et al. 2003; see also 
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Bergin et al. 2004). These authors speculate that the optically thin region arises from the 
presence of a companion that has carved out a gap in the disk (see also Marsh & Mahoney 
1992). 

In comparison, the [Nell] emission component centered on the systemic velocity in GM 
Aur has a HWZI of ~ lOkms -1 . If the width of the emission arises from Keplerian rotation, 
the emission arises from > 4 AU, given the stellar mass and inclination of GM Aur (Dutrey 
et al. 1998; Simon et al. 2000). This component may therefore arise from within the optically 
thin region (5-24 AU) in the interpretation of the SED favored by Calvet et al. (2005; see also 
Dutrey et al. 2008), or from the inner region of the outer disk, if the outer disk extends in to 
the smaller radii favored by Rice et al. (2003). Figure 5 shows a model profile constructed 
using the P(v K ) distribution of Meijerink et al. (2008) where the emission extends in to 4 AU 
(long dashed green line). Higher signal-to- noise ratio data are needed to determine if this 
model profile provides a better fit than one in which the inner radius extends in to much 
smaller radii (e.g., 0.25 AU; short dashed blue line). 

The redward component of the [Nell] emission extends to ~ 33 — 55kms _1 from the 
stellar velocity. If the [Nell] line profile is symmetric across the stellar velocity and arises 
from Keplerian rotation, the redward emission component can be inferred to arise from disk 
radii 0.16—0.5 AU. This component may represent emission from a gaseous inner disk close to 
the star. Such a gaseous disk is expected to be present since GM Aur experiences continued 
stellar accretion at a rate typical of T Tauri stars (~ lO~ 8 M yr _1 ; Najita et al. 2007). 

5. Summary and Future Directions 

The observations reported here bring to a grand total of four the number of T Tauri 
[Nell] emission sources that have been studied at high spectral resolution. In the high 
accretion rate system T Tau, the [Nell] emission is associated primarily with outflow activity 
(van Boekel et al. 2009). In the three lower accretion rate cases studied thus far (TW Hya, AA 
Tau, GM Aur), the [Nell] emission appears to be centered at the stellar velocity, consistent 
with a disk origin for the emission. The [Nell] emission is broader in the non-face-on systems 
(AA Tau, GM Aur) than in the face-on TW Hya, a result that is also consistent with a disk 
origin for the emission. 

In the non-face-on systems, the [Nell] emission width is narrower than (or comparable 
to) the CO emission width. If the widths of both diagnostics are dominated by Keplerian 
rotation, this suggests that the [Nell] emission arises from larger disk radii on average than 
the CO emission. This is consistent, in general, with models of the ionization and thermal 
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structure of T Tauri disks (e.g., Glassgold et al. 2007; Meijerink et al. 2008). 

The observed [Nell] line profiles show little evidence for an origin in a disk photoevap- 
orative flow. This is particularly so in the case of GM Aur, where the central component 
of the [Nell] emission is narrow enough to plausibly arise in a photoevaporative flow, but 
the profile is centered at the stellar velocity rather than showing the 5 — lOkms -1 blueshift 
that is predicted for emission arising in a photoevaporative flow (Alexander 2008). Pho- 
toevaporative flows may be present in these systems but contribute a small fraction of the 
[Nell] emission. More recent work by Pascucci & Sterzik (2009) finds stronger evidence for 
photoevaporative flows in other [NeII]-emitting systems. 

The equivalent width of the [Nell] emission we detect is less than that of the spectrally 
unresolved [Nell] feature in the Spitzer spectra of the same sources. Variability in the [Nell] 
emission or the mid- infrared continuum, a spatially extended [Nell] component, or a very 
(spectrally) broad [Nell] component might account for the difference in the equivalent widths. 
Further work is needed to understand the origin of this discrepancy. 

These results illustrate the ability of high resolution spectroscopy to probe the origin of 
the [Nell] emission from T Tauri stars. Further measurements of resolved [Nell] line profiles, 
for both classical T Tauri stars and transition objects, are needed to determine whether the 
results obtained here apply to the majority of [NeII]-emitting T Tauri stars. In addition, 
higher sensitivity line profiles than those reported here would be useful to probe the origin 
of the [Nell] emission in detail. 
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Fig. 1. — Spectrum of [Nell] emission from AA Tau (histograms) plotted as a function of 
heliocentric velocity and smoothed by 3 pixels. The vertical scale is flux density in units of 
Jy. Spectra obtained on 29 Oct 2007 and 30 Oct 2007 (top and middle panels) are shown 
along with the combined spectrum (bottom panel). The light lines in the upper two panels 
indicate the noise per pixel calculated by the TEXES pipeline at each point in the spectrum. 
The gaps in the spectra in the upper two panels reflect the gaps between the TEXES orders 
at this wavelength. The vertical dashed line in the bottom panel indicates the stellar velocity. 
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Fig. 2. — Spectrum of [Nell] emission from GM Aur plotted as a function of heliocentric 
velocity and smoothed by 3 pixels. The vertical scale is flux density in units of Jy. The verti- 
cal dashed line indicates the stellar velocity, and the horizontal line indicates the continuum 
level. The light lines indicate the noise per pixel calculated by the TEXES pipeline at each 
point in the spectrum. In addition to the [Nell] emission component centered on the stellar 
velocity, another emission component appears to be present at +40kms _1 from the stellar 
velocity. 
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Fig. 3. — Top: CO emission profile from AA Tau in the 4.9 /im region, normalized to the 
continuum. The vertical dashed line indicates the stellar velocity (i>heiio = 16.5 km s" 1 ). The 
line profile shown is the average of the v =1-0 P25 through P28 lines. Regions of poor telluric 
correction were excluded from the average. Bottom: [Nell] line profile of AA Tau from the 
combined spectrum shown in Figure 1, normalized to the continuum. 
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Fig. 4. — Top: CO emission profile from GM Aur in the 4.7 ^tm region, normalized to the 
continuum. The vertical dashed line indicates the stellar velocity (fViio — 14.7 km s -1 ). The 
line profile shown is the average of the u=l— P8 through P12 lines. In constructing the 
profile, the stellar photospheric contribution to the spectrum was first removed and regions 
of poor telluric correction were excluded from the average. Bottom: [Nell] line profile of GM 
Aur from the spectrum shown in Figure 2, normalized to the continuum. 
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Fig. 5. — [Nell] emission profile of GM Aur, compared with synthetic [Nell] disk emission 
profiles from an X-ray irradiated disk (Meijerink et al. 2008; see text for details). The model 
profiles assume Keplerian rotation and the stellar mass (O.85M ) and inclination (i = 54 
degree) of GM Aur reported by Simon et al. (2000). The emission extends from an outer 
radius of ~ 25 AU, beyond which the [Nell] emission is insignificant, inward to 0.25 AU 
(short dashed blue line) or 4 AU (long dashed green line). A vertical scaling factor is applied 
to each model to facilitate comparison with the observed profile. 
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Table 1. [Nell], H 2 , and CO Line FWHM 



Source 


i 


a(Nell) 


<r(H 2 ) 


<r(CO) 


Refs. 






kms -1 


km s _1 


kms -1 




AA Tau 


75 


70 




145 


B99, Here 


GM Aur 


54 


14 




30 


S00, S07, Here 


TW Hya 


4 


21 


14 


8 


P08, H07, S07 


BP Tau 


< 50 




57 


70 


A02, N03 


DF Tau 


80 




27 


65 


A02, N03 








23 




H06 


RW Aur 


40 




52 


80-250 


A02, N03 


V836 Tau 


65 




24 


130 


H06, N08 



a A02=Ardila et al. 2002; B99=Bouvier et al. 1999; 
H06=Herczeg et al. 2006; H07=Herczeg et al. 2007; N03=Najita 
et al. 2003; N08=Najita et al. 2008; P08=Pontoppidan et al. 2008; 
S00=Simon et al. 2000; S07=Salyk et al. 2007 



